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Abstract 

We report recent results of a WIMP dark matter search experiment using 310g of 
CaF2(Eu) scintillator at Kamioka Observatory. We chose a highly radio-pure crystal, 
PMTs and radiation shields, so that the background rate decreased considerably. We 
derived limits on the spin dependent WIMP-proton and WIMP-neutron coupling 
coefficients, and a^. The limits excluded a part of the parameter space allowed 
by the annual modulation observation of the DAM A Nal experiment. 
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1 Introduction 



There is substantial evidence that most of the matter in our Galaxy must 
be dark matter that exists in the form of Weakly Interacting Massive Par- 
ticles (WIMPs) [1]. WIMPs arc thought to be non-baryonic particles, and 
the most plausible candidates for them are the lightest supersymmetric par- 
ticles. WIMPs can be directly detected through elastic scattering with nuclei 
in radiation detectors. Using various detectors, many groups have performed 
experiments for the detection of WIMPs. 

For the direct WIMP detection, two kinds of interactions need to be consid- 
ered: the axial-vector (Spin-Dependent, SD) interaction and the scalar (Spin- 
Independent, SI) interaction. WIMPs couple to the spin of the target nuclei in 
the SD interaction while they coherently couple to almost all nucleons of the 
target nuclei in the SI interaction. For the SD interaction, ^^F is one of the 
most favorable nuclei to detect WIMPs because of its large nuclear spin and 
100% natural abundance. Furthermore, the spin expectation values of protons 
and neutrons in ^^F have the opposite signs while those values in other nuclei 
usually have the same signs [2-6]. Therefore, experiments with a ^^F-target 
can set complementary limits to those with nuclei whose spins have the same 
sign. 

Several direct WIMP searches using ^^F-based detectors, such as bolometers 
[7,8], scintillators [9,10] and superheated droplet detectors (SDDs) [11,12], 
have already been performed. Our group had also carried out WIMP search 
experiments using LiF and NaF bolometers [7,8]. Although these results set 
complementary limits to those of Nal(Tl) experiments, the background rates 
below 20 keV could not be decreased sufficiently. The sources of the main 
background were thought not to be any radiation but to be intrinsic and in- 
strumental noises in the bolometers. As an alternative method, we have carried 
out experiments using CaF2(Eu) scintillators. CaF2(Eu) is known as a useful 
scintillator for WIMP searches because of its high light yield (19000 pho- 
tons/ Me V) and has already been used in several experiments. However, these 
experiments show rather high background rates (^10 counts/k.e.e./day/kg ^ ) 
in comparison with the Nal(Tl) experiments (~ 1 counts/k.e.e./day/kg). In 
addition, its light yield is about 50 % of Nal(Tl) , so that the energy thresh- 
old of CaF2(Eu) detectors is higher than Nal(Tl) detectors. Therefore, these 
CaF2(Eu) experiments set less stringent limits to the SD interaction than the 
Nal(Tl) experiments. In order to reduce the background rates and to improve 
the energy threshold of our experiments, we have tried to eliminate radioactive 



^ In this letter, we use the unit k.c.c. or kcV electron equivalent for nuclear recoil 
energy since the scintillation efficiency for nuclear recoils is less than electron recoils 
due to high ionization loss. 
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sources of 7-rays and studied detection efficiency near tlie energy tliresfiold. 

fn this letter, the results of our new experiment with a CaF2(Eu) scintillator 
are presented. 



2 Experimental setup 

The detector consisted of a CaF2(Eu) crystal with a mass of 310g. It was 
installed in Kamioka Observatory (2700 m.w.e.). To reduce the background, 
low radioactive PMTs (Hamamatsu R8778) and a low radioactive radiation 
shield were used. The experimental setup is shown in Fig 1. The crystal was 
produced from the same stock of CaF2 raw powder as used by the CANDLES 
experiment [14] and low radioactive EuFa powder. The low radioactive PMT 
was developed by the XMASS experiment [13] and thought to be the least 
radioactive among available PMTs. Two PMTs were attached to the crystal 
through 5cm-long quartz light guides. The photoelectron yield of the scin- 
tillator was about 4 photoclcctrons/k.c.c. at 60 k.e.e. The radiation shield 
consisted of 5cm of highly pure copper (99.9999%), 10 cm of OFHC copper, 
15 cm of lead, and 20 cm of polyethylene. The highly pure copper was supplied 
by Mitsubishi Materials. The whole setup was separated from the mine air by 
two layers of EVOH sheet [15] so that radon gas in the air could not come into 
the detector. The radon free air generated in Kamioka observatory was sent 
into the detector as shown in the figure. In addition, the radon concentration 
around the detector was continuously monitored by a radon detector. It was 
kept at about 30 mBq/m^ throughout the experiment. 

For the data acquisition, a dual trace digital oscilloscope was used to record 
wave forms of PMT pulses. The pulses of each PMT are directly sent to the 
oscilloscope without any amplifier then daisy-chained to a trigger system with 
a tee connector. The impedance of each oscilloscope input was set to IMfl and 
the cables are terminated at the inputs of the trigger system. The wave forms 
were digitized at the oscilloscope at a rate of 1 GS/s for a total digitization 
time of 10 fj,s and sent to a PC for an off-hne analysis. 

In the trigger system, the pulses were amplified by PMT amplifiers and sent to 
low threshold discriminators. Thresholds of the discriminators were set to 1/4 
of the mean pulse hight of single photoelectron events. The single photoelec- 
tron events were obtained by irradiating the PMT with continuous feeble fight 
from an LED. The detection efficiency of each PMT with the discriminator 
threshold for a single photoelectron event was obtained to be 0.8 by comparing 
the discriminator output counts and the total number of single photoelectron 
events; the latter was estimated by assuming a normal distribution for the 
charge distribution of the single photoelectron events. Parameters of the nor- 
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mal distributions were determined in a special measurements by recording the 
events with low-threshold internal triggers of the oscilloscope to include the 
low energy tail of the spectrum. The charge distribution of the internally trig- 
gered events was well fitted by the normal distribution. Single counting rates 
of PMTs were typically 20 Hz. 

The discriminator signals were then sent to the coincidence circuit. It gener- 
ated a trigger pulse for the oscilloscope when both of the PMTs gave signals 
in coincidence within l//s. It reduced the number of background events due to 
dark noise of the PMTs. 

The trigger efficiency of the data acquisition system was calculated as a func- 
tion of a number of photoelectrons by the Monte-Carlo analysis. In this anal- 
ysis, we generated a certain number of photoelectrons in a time sequence 
according to the scintillation decay time of CaF2(Eu). We assumed each pho- 
toelectron reached in each PMT evenly. The detection probability of a pho- 
toelectron was set to 0.8. Then, we measured an interval between arbitrary 
photoelectrons of each PMT in the event. The fraction of events that gener- 
ated at least one output signal in both of the discriminators within 1 /is was 
taken to be the trigger efficiency. At 2 k.c.c., 8 photoelectrons were generated 
and the trigger efficiency was calculated to be 0.93. Threshold energy for the 
following analysis was set at 2 k.e.e. 



3 Measured spectra 

The WIMP observation was carried out from March to May 2005. Measured 
spectra are shown in Fig 5. The energy scale of the spectra is defined to be 
proportional to the number of photoelectrons. In the energy region below 10 
k.e.e., background events due to Cherenkov photons dominate the event rates. 
The Cherenkov events were produced in the light guides by Compton electrons 
caused by background 7-rays. We used a pulse shape discrimination (PSD) 
technique to eliminate these events. Cherenkov photons are observed as fast 
pulses (< 10 ns) while scintillation photons as slow pulses (~ 1 yus). Therefore, 
we used the ratio of the integral charge of the partial pulse shape period (0-30 
ns) to that of the total pulse shape period (0-10 ^s) as a discrimination pa- 
rameter in the off-line analysis. To set the discrimination cut, we compared low 
energy Compton events which cannot produce Cherenkov photons (obtained 
using 122 keV 7-rays from ^'^Co) and Cherenkov events (obtained using 1133 
and 1333 keV 7-rays from ^°Co). In addition, we eliminate events in which 
observed charge of the two PMTs shows high asymmetry because they are 
thought to be the remaining Cherenkov events or electric noise events. To es- 
timate the efficiency of the off-line event selection, the low energy Compton 
events were used. We calculated a fraction of the events which remained after 
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the selection as a function of energy. It was 0.74 at 2 k.e.e. Consequently, the 
event selection reduced the count rates to less than 10 counts/k.e.e./day/kg 
in the energy region between 2 and 10 k.e.e. 



4 Limits on the WIMP-nucleon interaction 



Using the measured spectrum, limits on the SD WIMP-nucleon interaction 
were derived. The limits were calculated using the same manner as described 
in Ref. [16]. To derive the upper limits, we conservatively assumed all events 
to be nuclear recoils caused by the WIMP scattering. The astrophysical and 
nuclear parameters used to obtain the limits are listed in Table 1. Only the 
contribution of ^^F was considered while ^°Ca had negligible contribution to 
the SD interaction because ^°Ca consisted of even protons and even neutrons. 
The scintillation efficiency /q of the ^^F recoil in a CaF2(Eu) was measured 
by other experiments. Although all of these experiments show that the scintil- 
lation efficiency tends to increase in lower energy region, we choose conserva- 
tively a constant value /q = 0.11. As was done in Ref. [17], ^^F recoil energy is 
calculated with /q and the interpolated energy scale between and 60 k.e.e.: 

= 7 7 (Ij 

-t^eOkeV Jq 



where E'n is the ^^F recoil energy and L is the scintillation yield of an event. 
For the SD interaction, the WIMP-nucleus cross section is written as 

SD 32G|//^_N/ /c \ I /c \^2'^+ -'- /o\ 

where Gp is the Fermi coupling constant, /U^-n is the WIMP-nucleus reduced 
mass, (iS'p(N)) and (<S'n(N)) are the expectation values of the proton and neutron 
spins within the nucleus and J is the total nuclear spin [18]. 

The upper limit of the WIMP-nucleus scattering cross section cry2j^_N is ob- 
tained by the experiment. From Eq. (2), limits in the ap-On plane is written 
with afP^_^: 



J 32GpyU^_f^ 



As in Ref. [19], limits on a single nucleon interaction are generally calculated: 
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where {Sp) = (5'„) = i, o-f°^_p(N) and crj^^^.^^N) are the WIMP-proton 
and WIMP- neutron scattering cross section hmits when an(/S'n(N)) — and 
ap('Sp(N)> = in Eq. (3), respectively. Figure 3 shows (T,f°^_p(N) and (7i^^^_^(n) 
derived from our experiment as a function of WIMP mass M^. Results of other 
experiments are also shown in the figure. 

The hmits in the Op-On plane for WIMP mass = 50 and 200 GeV are 
shown in Fig 4. The limits of other experiments are also shown. The outside 
of the two solid lines is excluded by our experiment. The region between the 
two ellipses is allowed by the annual modulation observation of the DAMA 
experiment [24]. Our results exclude a part of the DAMA allowed region by 
means of odd-proton target. Experiments with odd- neutron targets such as the 
CDMS [21] also set similar limits. Our results are comparable to the recent 
results of PICASSO [12] experiment which uses SDDs of fluorocarbon, C4F10. 



5 Conclusion 

The WIMP search experiment using 310g of CaF2(Eu) was carried out. The 
highly pure copper shield and the quartz light guides were used to eliminate 
7-rays from the outer materials. Pulse shape discrimination effectively elim- 
inated Cherenkov events, which dominated the count rates below 10 k.e.e. 
As a result, count rates were lower than 10 counts/k.e.e. /day/kg between 2 
and 10 k.e.e. We obtained the limits on the WlMP-nucleon spin dependent 
interaction in terms of the WlMP-nucleon coupling coefficients ttp and an. 
Our results excluded a part of the parameter region allowed by the annual 
modulation measurement by the DAMA Nal experiment. 
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Astrophysical parameters 
Dark matter density 
Velocity distribution 
Velocity dispersion 
Solar system velocity 



0.3 GeV/cm^ 
Maxwellian 
220 km/s 
232 km/s 



Nuclear parameters of -"^^F 
Total spin 

Spin expectation value of protons 
Spin expectation value of neutrons 



1 

2 

0.441 
-0.109 



Table 1 

Astrophysical and nuclear parameters used to obtain the limits. 
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Fig. 1. Schematic view of the experimental setup. 
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Fig. 2. Measured spectra before and after the event selection is applied. The energy- 
scale is normalized to a calibration point of 60 keV electron recoil. 
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Fig. 3. Recent limits on the SD WIMP-proton cross section and WIMP-neutron 
cross section. Results of NAIAD [20], PICASSO [12], CDMS[21], DAMA Nal[22], 
DAMA Xe[23] and our experiments are shown. 
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Fig. 4. Limits in the Op-On plane for = 50 GeV (upper) and 200 GeV (lower). The 
region between two solid lines is allowed by this experiment. Results of NAIAD [20], 
PICASS0[12], CDMS[21], DAMA Nal[22] and DAMA Xe[23] experiments are also 
shown. 
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